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I.  Introduction 

This  is  the  Final  Technical  Report  on  work  performed  under  the 
support  of  the  Air  Force  Office  of  Scientific  Research  under  Grant  AFOSR 
88-0027  for  the  period  October  1,  1989,  to  September  30,  1990.  The 
objective  of  this  work  was  to  carry  out  basic  physics  research  on  the 
generation  of  coherent,  tunable  radiation  from  relativistic  electron  beams 
and  exploring  means  of  improving  performances. 

II.  Summary  of  Work  Accomplished 

This  work  covered  computer  simulations  and  theoretical  studies  of 
radiation  generation  in  three  different  free  electron  microwave  devices:  the 
cyclotron  autoresonance  maser  (CARM),  the  backward  wave  oscillator 
(BWO),  and  the  magnicon. 

In  investigating  the  effects  of  magnetic  field  on  the  output  power  of  a 
plasma-filled  Backward  Wave  Oscillator,  we  found  that  within  a  certain  range 
of  magnetic  field  the  growth  rate  of  beam-plasma  cyclotron  interaction  is 
significantly  larger  than  the  conventional  backward  wave  oscillation.  This 
slow  plasma  cyclotron  wave  is  emitted  from  plasmas  only  if  it  also  satisfies 
the  dispersion  relation  of  a  ripple  wall  waveguide  which  couples  the  slow 
wave  to  fast  wave.  This  observation  can  be  employed  to  more  correctly 
interpret  previous  experiments  carried  out  in  the  Soviet  Unions  which 
displayed  a  strong  dependence  of  output  power  on  the  magnetic  field. 

We  have  carried  out  computer  simulations  of  a  100  GHz  electron 
cyclotron  autoresonance  maser  oscillator  that  was  designed  for  an 
experiment  at  the  Naval  Research  Laboratory.  The  results  indicate  an 
intrinsically  low  efficiency  for  the  experimental  parameters.  Enhancement 
of  the  efficiency  by  varying  the  beam  current,  interaction  length,  magnetic 
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field, and  tapering  the  magnetic  field  are  investigated.  It  is  demonstrated 
that  an  efficiency  above  20%  may  be  obtained  by  reducing  or  tapering  the 
magnetic  field. 

In  simulating  the  performance  of  a  Magnicon,  we  found  that  the 
saturation  mechanism  for  the  uniform  magnetic  field  case  is  due  to  the 
phase  detuning  between  the  electron  beam  and  electromagnetic  wave.  On 
the  other  hand,  the  conversion  of  electron  kinetic  energy  to  radiation  in  the 
tapered  magnetic  field  case  is  terminated  only  if  the  electrons  deplete  their 
entire  transverse  energy. 

A  unifed  theory  is  developed  for  free  electron  lasers,  cyclotron 
autoresonance  masers,  and  ion  channel  lasers.  Based  on  this  theory, 
parametric  studies  have  been  carried  out  to  determine  when  different 
bunching  mechanisms  enhance  or  interfere  with  each  other. 

Below,  we  summarize  the  scope  and  main  findings  of  these 
investigations.  More  detailed  information  are  contained  in  the  appended 
reports. 

II.  1  Emission  of  Plasma  Cyclotron  Waves  in  Plasma-Filled  BWO 

(Appendix  1) 

Previously,  there  were  two  experiments  in  relativistic  plasma-filled 
BWO  which  studied  the  dependence  of  output  power  on  the  magnetic  field. 
Both  results  illustrated  that  the  output  power  peaked  at  a  certain  magnetic 
field  and  they  all  conjectured  that  the  resonance  was  due  to  nf2b  =  k0V0. 
However,  the  experimental  parameters  used  did  not  satisfy  this  resonance 
condition.  Both  experiments  employed  a  small  solid  beam. 

We  have  carried  out  computer  simulations  using  the  experimental 
parameters  to  elucidate  the  physics  mechanism  which  gives  rise  to  the 
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resonance  behavior.  A  vacuum  BWO  simulations  with  Q5  =  10coc  have  been 
carried  out.  No  instability  (conventional  BWO)  has  been  observed  for  a  solid 
beam  with  Rb  =  0.3Rq  which  was  used  in  the  early  experiment.  In  the 
presence  of  the  dense  background  plasma  (cop  =  0.45coc)  and  using  Qp=  1.5 
7coc,  the  slow  plasma  cyclotron  wave  is  observed  to  be  unstable.  This  wave  is 
mainly  supported  by  the  background  plasma  and  it  can  only  leak  out 
evanescently  if  no  wave  conversion  mechanism  is  provided.  The  detector 
placed  outside  the  plasma  will  be  very  difficult  to  observe  this  instability.  By 
varying  the  magnetic  field,  the  intersection  between  the  beam  slow  space 
charge  mode  and  the  plasma  cyclotron  mode  can  be  brought  to  coincide 
with  the  dispersion  curve  of  the  electromagnetic  branch  in  a  ripple 
waveguide.  The  excitation  of  the  fast  wave  branch  of  the  TM  waveguide 
mode  is  observed  only  if  fip  is  in  the  neighborhood  of  1.32coc.  This  result 
reproduces  quite  well  the  experimental  curve  which  displays  the  output 
power  versus  magnetic  field.  Ihe  observed  mechanism  will  have  important 
implications  in  building  a  medium  power  beam-plasma  microwave  device 
which  will  be  discussed  in  section  III.  1 . 

II. 2  Cyclotron  Autoresonance  Maser  Oscillator 

(Appendix  2) 

In  a  CARM  amplifier  an  injection  source  is  required  which  in  some 
frequency  region  may  not  be  readily  available.  Recently,  attempts  have  been 
made  at  the  Naval  Reserach  Laboratory  to  develop  a  100GHz,  24MW  CARM 
ocillator  operating  with  20%  efficiency.  Initial  results  show  disappointingly 
low  efficiency  with  mode  competition  suspected  to  be  the  cause. 

The  main  findings  for  our  simulation  results  are:  Even  in  the  absence 
of  mode  competition  (priming)  the  efficiency  at  the  experimental 
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parameters  is  very  low,  only  about  3%.  Higher  efficiency  may  be  achieved  by 
either  using  a  lower  magnetic  field  or  tapering  the  magnetic  field.  By 
tapering  the  magnetic  field  simulations  using  a  higher  beam  current  of  400 
amps  show  it  is  possible  to  attain  an  efficiency  above  20%  for  an  electron 
beam  with  no  velocity  spread,  and  above  10%  with  the  experimental 
anticipated  maximum  velocity  spread  of  3%. 

No  gyrotron  mode  competition  was  observed  in  any  of  the  simulations. 
The  efficiencies  were  thus,  obtained  under  essentially  ideal  single-mode 
operation.  The  absence  of  gyrotron  mode  competition  could  be  due  to  our 
priming  the  oscillator  at  the  CARM  frequency,  which  allows  the  CARM  mode 
to  start  at  a  much  higher  power  relative  to  the  gyrotron  mode  (which  has  to 
start  from  the  intitial  noise).  We  find  that  without  the  priming,  the  initial 
broad  band  kz-  spectrum  of  the  rf  field  from  noise  will  evolve  into  two  peaks 
centered  at  the  CARM  and  gyrotron  resonances,  and  the  rf  field  takes  a  long 
time  to  grow.  This  confirms  the  possibility  that  priming  could  be  beneficial 
in  discriminating  against  unwanted  longitudinal  as  well  as  transverse  modes. 

II. 3  Magnicon 

Most  microwave  free  electron  devices  start  with  random  phases 
between  electrons  and  electromagnetic  wave  and  rely  on  electron  bunching 
to  efficiently  generate  microwave  radiation.  On  the  other  hand,  there  are 
generation  mechanisms  which  require  no  electron  bunching  by  injecting  the 
electrons  into  the  cavity  with  optimal  single  phase,  such  as  the  Magnicon 
which  has  recently  been  shown  to  produce  2.6MW  output  microwave  pulse 
of  30ps  at  0.915GHz  with  85%  efficiency.  The  output  cavity  works  in  the 
TM  no  mode  and  the  device  operates  on  gyroresonance  (co0  =  O)  without  the 
Doppler  shift  term,  and  the  electron  entrance  into  the  output  resonator  is 
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varied  continuously  in  synchronism  with  the  wave  rotation.  The  energy 
conversion  from  electrons  to  the  wave  is  determined  from  an  approximate 
balance  relation  of  the  Lorentz  forces  between  eEM=  -  ecorBxand  evxBi  = 
ewrB±.  Thus,  during  the  wave-particle  interaction,  the  electron  transverse 
energy  is  converted  into  longitudinal  energy  and  dissipated  into  En.  The 
axial  momentum  should  remain  constant  which  implies  that  as  the  electron 
loses  energy,  its  axial  velocity  increases. 

We  have  developed  a  self-consistent  large  signal  code  to  study  the 
Magnicon  interaction.  The  preliminary  simulation  results  are  shown  in  Fig. 
1  and  2  using  the  TMno  mode,  and  the  following  parameters:  a  =  v1/vn=  4, 

Rg/PL=  1.  300kV  and  60  A  beam  (off-axis  encircling  beam).  The  interaction 
is  optimized  by  varying  the  axial  interaction  length  and  tapering  the  axial 
magnetic  field  (B0).  The  best  result  we  attain  so  far  uses  VB0/B0  =  0.2/L  (L  = 
44  cm).  The  time  evolution  of  the  electron  transverse  velocity  in  a  wave 
field  near  the  device  axis  and  a  resonant  axial  uniform  magnetic  field  can  be 
expressed  as 


V±  = 


v±0  +  ^HF  VZ0 


t2  -  2QHFVzot  sin<|) 


1/2 


(1) 


where  O  =  k  -  <j>,  and  k  is  the  electron  entrance  phase  and  <>  is  the  rf  phase. 
Equation  (1)  shows  that  when  O  =  7t/2,vx  decreases  linearly  in  time.  Figure  1 
display  the  time  evolution  of  the  electron  energy  (yR)  when  it  comes  out  of 

the  cavity  and  the  cavity  field  energy  (Wf)  which  demonstrates  that  tapering 
the  magnetic  field  increases  the  efficiency  from  65%  to  85%.  Figure  2 
shows  the  spatial  evolution  of  the  test  electron  energy  and  its  phase  relative 
to  the  wave  (0).  In  the  simulation  the  electromagnetic  wave  grows  from  the 
background  noise.  For  the  untapered  case,  0  starts  at  around  k/2  which  is 
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Comparison  between  tapered  and  untapered  Magnicon.  Time  evolution  of 
(a)  electron  energy  at  the  exit  plane,  (b)  resonator  field  energy. 
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the  optional  entrance  phase  (Eq.  1)  and  becomes  about  n  when  the 

— ^  — > 

conversion  process  terminates  (V^x  B±  =  0).  The  change  in  0  is  primarily 
due  to  the  increase  in  the  electron  relativistic  cyclotron  frequency  as  it  loses 
energy.  This  detuning  can  be  compensated  by  tapering  the  magnetic  field. 
In  this  case  the  conversion  process  ends  when  the  electron  transverse 
energy  vanishes.  Although  the  electron  starts  with  93%  of  its  energy 
residing  along  the  transverse  direction,  the  optimized  efficiency  reaches 
only  85%.  This  is  because  during  the  conversion  process  P|i  =  m0Yvn 
remains  constant.  As  y  decreases  Vn  has  to  increase.  We  continue 
investigating  the  Magnicon  and  the  investigation  will  be  extended  to  include 
the  harmonic  Magicon. 

II. 4  Unified  Theory  for  CARM.  ICL.  and  FEL 

(Appendix  3) 

Recently,  a  new  radiation  source,  the  Ion-Channel  Laser  (ICL)  was 
conceived  at  UCLA.  In  this  device,  a  relativistic  electron  beam  is  injected 
into  an  underdense  plasma  (cop  <  <%)•  The  beam  front  pushes  out  the  plasma 
electrons  leaving  an  ion  channel.  The  ion  focusing  force  causes  the  beam 
electrons  to  oscillate  about  the  axis  and  plays  a  similar  role  as  the  magnetic 
field  in  a  CARM.  Radiation  can  be  generated  ranged  from  microwaves  to 
soft-X-rays. 

We  have  developed  a  unified  theory  to  describe  the  amplification 
mechanism  for  CARM,  ICL,  and  FEL.  The  theory  assumes  that  the 
dependence  of  electron  gyrofrequency  on  its  beam  energy  is  in  the  form  of 
yQ  where  q  =  1,  and  0,  respectively  for  CARM,  ICL.  and  FEL.  The  main 
findings  are  :  Azimuthal  bunching  dominates  over  axial  bunching  if  q(l- 


Voz/Vph)  >  ("T"  -  ~)  and  vice  versa.  If  voz  >  vph  or  c2  <  voz  vph,  the  two 
vPh  Ph 

bunchings  reinforce  each  other  rather  than  offsetting  which  gives  rise  to  no 

auto  resonance.  In  order  to  have  wave  amplification,  the  sign  of  mismatch 

has  to  be  the  same  as  the  sign  of  bunching  parameter.  For  CARMs  and  ICLs 

the  relativistic  mass  variation  effect  on  the  electron  axial  velocity  may 

overwhelm  the  "v  x  B  force  so  that  the  chanr  in  axial  veloicty  is  opposite  to 

the  electron  energy  change  and  the  axial  bunching  term  becomes  negative. 

For  FELs  the  bunching  arising  from  the  relativistic  mass  variation  is  smaller 

than  that  of  the  ponderomotive  force  induced  by  the  wiggler  and  the  wave  so 

the  FEL  bunching  paramerer  is  always  negative.  All  these  devices  can 

2 

upshift  the  electron  undulation  frequency  by  a  factor  of  2yoz. 
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It  is  found  through  computer  simulations  of  plasma-filled  backward-wave  oscillators  that  within  a  cer¬ 
tain  range  of  magnetic  field  the  growth  rate  of  beam-plasma  cyclotron  interaction  is  significantly  larger 
than  the  conventional  backward-wave  oscillation.  This  slow  plasma  cyclotron  wave  is  emitted  from  plas¬ 
mas  only  if  it  also  satisfies  the  dispersion  relation  of  a  rippled-wall  waveguide  which  couples  the  slow 
wave  to  fast  wave.  This  observation  may  be  employed  to  interpret  previous  experiments  which  displayed 
a  strong  dependence  of  output  power  on  the  magnetic  field. 


PACS  numbers:  52.60. -Hi,  52.65. +z 

Recently,  a  considerable  amount  of  effort1'5  has  been 
undertaken  towards  the  production  of  ultrahigh -power 
microwave  sources  based  on  backward-wave  oscillator 
(BWO)  configurations,  in  which  a  rippled-wall  slow- 
wave  structure  provides  a  very  efficient  way  of  converting 
the  kinetic  energy  of  relativistic  electron  beams  into  mi¬ 
crowave  radiations. 

In  order  to  raise  the  relativistic  beam  current  and 
hence  the  output  power  of  relativistic  electronic  devices, 
a  background  plasma  is  sometimes  employed  to  provide 
space-charge  neutralization.  The  plasma  may  also  in¬ 
duce  some  other  effects  which  could  influence  the  out¬ 
come  of  backward-wave  oscillators.  Experimentally,  it 
has  been  observed6  that  the  efficiency  can  be  increased 
by  a  factor  of  8  when  the  plasma  frequency  (Op^Q.lcko 
(ko  is  the  rippled-wall  wave  number  and  c  is  the  speed  of 
light)  and  the  enhancement  mechanism  was  attributed  to 
the  excitation  of  beam-plasma  backward-wave  instabili¬ 
ty.  Very  recently,  computer-simulation  results7  have 
demonstrated  that  a  more  plausible  explanation  for  the 
efficiency  enhancement  is  that  a  dense  background  plas¬ 
ma  tends  to  reduce  the  phase  velocity  of  the  most  unsta¬ 
ble  mode  in  a  BWO  and  causes  the  beam  to  convert 
more  energy  to  the  wave. 

The  investigations  described  above  all  impose  a  strong 
magnetic  field  to  confine  the  electron  motion  to  only  one 
dimension.  In  vacuum  BWO’s,8  the  output  powers  were 
observed  to  depend  on  the  magnetic  field.  The  beam  cy¬ 
clotron  resonance  absorbed  the  BWO  emission  and  intro¬ 
duced  a  dip  in  the  curve  of  the  output  power  versus  the 
magnetic  field.  This  issue  will  not  be  treated  in  this  pa¬ 
per.  On  the  other  hand,  experiments9  in  relativistic 
plasma-filled  BWO  illustrated  that  the  output  power 
peaked  at  a  certain  magnetic  field  and  they  conjectured 
that  the  resonance  was  due  to  nftj,  “kok'o  (fit  is  the  rel¬ 
ativistic  beam  cyclotron  frequency,  n  is  the  harmonic 
number,  and  K0  is  the  beam  axial  velocity).  However, 
the  experimental  parameters  used  did  not  satisfy  this  res¬ 
onance  condition.  In  this  paper  we  will  demonstrate 
through  computer  simulations  that  the  excitation  of  slow 
plasma  cyclotron  waves  when  they  also  satisfy  the 
dispersion  relation  of  a  rippled-wall  waveguide  can  give 
rise  to  the  resonance  characteristics  observed  in  the  ex¬ 


periment.  This  mechanism  can  be  utilized  to  convert  the 
slow  wave  into  fast  wave,  and  thus  overcome  one  of  the 
major  obstacles  in  realizing  beam-plasma  interaction,  to 
build  a  compact  medium-power  microwave  device. 

The  dispersion  characteristics  of  a  plasma-filled 
rippled-wall  waveguide  [/?*— /?0(l+esin&oz),  e  is  the 
ripple  depth]  immersed  in  an  infinite  magnetic  field  have 
previously  been  shown6  and  confirmed  by  our  particle 
simulations.  They  consist  of  the  electromagnetic  branch 
of  the  rippled-wall  waveguide  and  a  plasma-wave 
branch.  The  effect  of  a  finite  magnetic  field  introduces  a 
plasma  cyclotron  wave  branch.  The  dispersion  relation 
of  the  plasma-wave  branches  in  a  plasma-filled  wave¬ 
guide  is  expressed  by10 

kc  2  _  \+a)p/(,(lp-a)2) 

°>c  1  —(Dp/a)1 

where  (oc  is  the  waveguide  cutoff  frequency,  o)p  is  the 
plasma  frequency,  and  Cip  is  the  plasma  cyclotron  fre¬ 
quency.  In  the  case  of  flp  >  tap,  the  uncoupled  disper¬ 
sion  relations  are  shown  in  Fig.  1  (a).  Notice  that  in  this 
situation,  the  frequency  of  the  plasma  cyclotron  wave 
branch  can  be  higher  than  the  cutoff  frequency  of  a 
plasma-filled  waveguide,  whereas  the  plasma-wave 
branch  is  always  lower.  A  finite  magnetic  field  will  also 
introduce  the  cyclotron  resonance  lines  associated  with 
the  relativistic  electron  beam  la>mmkVo±  ft Fig.  1(a)]. 
The  excitation  of  the  fast  cyclotron  wave11  (+ft(>) 
demands  a  certain  amount  of  beam  initial  transverse 
momentum,  whereas  the  slow  cyclotron  wave12  (  —  ft*) 
requires  a  slow-wave  structure. 

In  order  to  address  the  finite-magnetic-field  effects  on 
the  plasma-filled  BWO  output  power,  a  particle  wave¬ 
guide  code  whose  algorithms  have  been  described  previ¬ 
ously7  is  used.  Transversely,  only  the  TMoi  mode  is  re¬ 
tained  while  periodic  boundary  conditions  are  imposed 
along  the  axial  direction  so  that  the  nth  mode  wave  num¬ 
ber  is  defined  as  fc„“2nn/L,  where  L  is  the  system 
length  and  is  taken  to  be  256X,  (kg  is  the  spatial  grid 
size).  To  take  into  account  the  coupling  of  the  TEoi  and 
the  TMoi  modes  through  electron  gyration,  simulations 
which  included  both  modes  were  also  carried  out.  The 
results  only  differed  by  a  few  percent  from  the  single- 
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FIG.  I.  The  propagation  characteristics  of  a  plasma-filled 
rippled  waveguide:  (a)  the  dispersion  relation  with  a  finite 
magnetic  field.  (PW,  plasma  wave;  PCW,  plasma  cyclotron 
wave;  BSCW,  beam  space-charge  wave;  FBCW,  fast-beam  cy¬ 
clotron  wave.)  (b)  The  time  evolution  of  wave  energy  of  unsta¬ 
ble  modes  (n,  ”1.57<uf). 

mode  simulation.  To  facilitate  the  comparison  between 
experimental  and  simulation  results,  the  experimental 
parameters8  were  used  in  our  simulation.  The  beam  en¬ 
ergy  and  current  are  respectively  700  kV  and  4  kA.  The 
rippled-wall  waveguide  dimensions  are  the  following: 
/?o~13.9Xf,  Xx"0.1  cm,  €"0.1,  and  ko—k^  for  n*  16 
in  our  simulation  system.  Plasma  electrons  are  assumed 
to  fill  up  uniformly  the  waveguide  cross  section  up  to 
R-R0(\  ~e)  and  /?— 0.3/?o  for  beam  electrons.  All 
frequencies  are  normalized  to  coc.  The  plasma  density 
can  be  varied  and  it  will  affect  the  location  where  the 
output  power  is  peaked  with  respect  to  the  magnetic 
field.  The  following  simulation  results  are  for 
aip  —0.45<uf.  Simulations  using  Ro”22.5Xg  were  also 
carried  out.  The  resonance  behavior  of  the  output  power 
versus  the  magnetic  field  remains  the  same. 

Consider  a  solid  beam  without  any  initial  transverse 
momentum.  The  conventional  beam-forward  plasma- 
wave  instability  can  be  avoided  if  (op/(oc  <  yoYo/c  (yo  is 
the  beam  relativistic  factor),  which  is  satisfied  in  our 
simulation.  The  BWO  instability  (surface  wave)  is  most 
effective  when  a  hollow  beam  placed  very  close  to  the 
waveguide  wall  is  employed.  A  vacuum  BWO  simula¬ 
tion  with  ft*  “10a>c  was  carried  out.  A  very  weak  insta¬ 
bility  was  observed  for  a  solid  beam  with  R(,“0.3Ro 
which  was  used  in  the  early  experiment.9  No  significant 
growth  was  observed  over  the  time  of  the  simulation  run. 
As  was  demonstrated  in  our  early  work,7  using  a  large 


radius  beam,  the  growth  rate  estimated  from  the  simula¬ 
tion  result  agrees  quite  well  with  the  theoretical  predic¬ 
tion.  In  the  presence  of  the  dense  background  plasma 
and  using  ftp~1.57o>c,  the  time  evolutions  of  the  most 
unstable  wave  energy  are  shown  in  Fig.  1(b).  The 
growth  rate  estimated  from  simulation  results  is  <w, 
*"0.058o»f,  which  is  rather  strong  because  it  is  a 
volume-wave  interaction.  However,  this  slow  plasma  cy¬ 
clotron  wave  is  mainly  supported  by  the  background 
plasma  and  it  can  only  leak  out  evanescently  if  no  wave 
conversion  mechanism  is  provided.  With  the  detector 
placed  outside  the  plasma  it  would  be  very  difficult  to  ob¬ 
serve  this  instability.  The  frequency  spectrum  of  the  un¬ 
stable  mode  (n“13)  peaks  at  <u”1.62(uc.,  which  falls  on 
the  dispersion  curve  of  the  plasma  cyclotron  wave.  In 
this  case,  the  fast-electromagnetic-wave  component 
(n  ™  —  3)  is  relatively  small  in  comparison  with  the  plas¬ 
ma  cyclotron  wave. 

By  varying  the  magnetic  field,  the  intersection  between 
the  beam  slow  space-charge  mode  and  the  plasma  cyclo¬ 
tron  mode  can  be  brought  to  coincidence  with  the  disper¬ 
sion  curve  of  the  electromagnetic  branch  in  a  rippled- 
wall  waveguide.  In  conventional  BWO  interactions,  the 
resonance  conditions  {co“kYo  —  (Oblyi11  and  co ” [to7 
+  (*-ito)2c2]l/2l  are  usually  satisfied  simultaneously. 
These  expressions  indicate  that  the  slow  space-charge 
beam  mode  in  the  zeroth-order  Brillouin  zone  should  be 
coupled  with  the  TM  waveguide  mode  in  the  first-order 
Brillouin  zone  and  grow  together  with  the  same  growth 
rate.  To  illustrate  the  magnetic-field  resonance  effect,  a 
series  of  simulations  using  different  magnetic-field 
strengths  have  been  carried  out  and  the  results  arc 
displayed  in  Fig.  2.  The  excitation  of  the  fast-wave 
branch  of  the  TM  waveguide  mode  is  observed  only  if 
np  is  in  the  neighborhood  of  1.32<ac.  The  wave  energy 
time  evolutions  of  the  plasma  cyclotron  mode  (n  —  1 1 ) 
and  the  fast  electromagnetic-wave  mode  («"—  5)  are 
shown  in  Fig.  2(a).  Their  growth  rate  and  real  frequen¬ 
cy  are  respectively  <u,  — 0.065<uf  and  <z)r  “1.31 5a>c.  The 
peak  efficiency  of  the  interaction  has  achieved  about 
17%.  In  this  case,  the  fast-electromagnetic-wave  com¬ 
ponent  (n  —  —  5)  is  significantly  larger  than  the  plasma 
cyclotron  wave.  Since  the  plasma  electrons  participate 
actively  in  supporting  the  excitation  of  the  plasma  cyclo¬ 
tron  wave,  the  saturation  mechanism  will  involve  the 
trapping  of  a  portion  of  plasma  electrons  by  the  axial 
electric  field  of  the  excited  slow  wave.  This  has  been 
demonstrated  in  Figs.  2(b)  and  2(c),  which  plot  the  elec¬ 
tron  axial  velocity  versus  distance  at  mct  “0  and  at  the 
time  of  saturation  (a>fr“140),  which  clearly  illustrates 
the  active  participation  of  plasma  electrons  in  the  forma¬ 
tion  of  trapped  structures.  In  doing  so,  the  plasma  elec¬ 
trons  gain  significant  amounts  of  kinetic  energy  at  the 
expense  of  beam  energy  and  thus  reduce  the  overall 
efficiency  in  microwave  generation.  Figure  3(a)  shows 
the  time  evolution  of  the  beam  and  the  plasma  kinetic 
energy.  At  the  time  of  saturation,  the  beam  electrons 
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FIG.  2.  Simulation  results  of  a  plasma-filled  BWO  (np 
“1.32<uf):  (a)  The  time  evolution  of  wave  energy  of  unstable 
modes;  electron  axial  velocity  vs  distance  at  (b)  cod  "0  and  (c) 
cod  —  1 40. 


have  on  the  average  lost  30%  of  their  kinetic  energy. 
However,  almost  half  of  it  is  spent  in  increasing  the  plas¬ 
ma  kinetic  energy.  In  comparison  with  a  conventional 
BWO,  this  could  constitute  a  disadvantage  of  construct¬ 
ing  microwave  tubes  based  on  the  excitation  of  plasma 
waves.  The  wave  energy  versus  magnetic  field  is  plotted 
in  Fig.  3(b).  The  results  clearly  reveal  pronounced  peak 
at  nA,“1.32<u0.  Furthermore,  the  simulation  results  also 
show  that  for  the  absence  of  the  background  plasma  the 
peak  no  longer  exists.  In  the  resonance  region 
(np  —  1.25eoc  —  l.4<oc),  the  output  frequency  from  simu¬ 
lations  remains  relatively  constant  (<u,  “1.315coc).  This 
is  due  to  the  fact  that  the  plasma  cyclotron  wave- 
dispersion  curve  is  modified  when  the  unstable  mode  is 
also  a  normal  mode  of  the  rippled-wall  waveguide.  Out¬ 
side  the  resonance  region,  the  unstable  mode  frequency 
follows  the  plasma  cyclotron  wave-dispersion  relation. 
However,  in  these  cases,  the  wave  energy  is  mostly  dom¬ 
inated  by  the  slow  wave  and  saturated  at  low  levels. 
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FIG.  3.  Simulation  results;  (a)  The  time  evolution  of  beam 
and  plasma  energy  (np-1.32wr);  (b)  the  emitted  energy  vs 
magnetic  field. 

Based  on  the  fact  that  the  emission  can  be  detected  only 
if  a  fast  wave  is  generated,  the  experimentally  observed 
radiation  outside  the  resonance  region9  was  dominated 
by  the  weakly  growing  vacuum  BWO  instability.  This 
can  explain  why  the  output  frequency  measured  in  the 
experiment  is  independent  of  the  magnetic  field. 

These  findings  may  be  used  to  interpret  previous  ex¬ 
perimental  results  which  showed  that  the  emission  was 
peaked  at  Z?o"“3.85  kG  and  f)x  ™0.3  (p±  ~v±/c).  They 
attributed  that  in  addition  to  the  conventional  BWO  res¬ 
onances,  the  fast-beam  cyclotron  resonance  is  also  sat¬ 
isfied  [«"(&  —  A:o)Fo+nfl<,].  They  claimed  that  if 
koVo“ntlb,  the  fast-beam  cyclotron  resonance  becomes 
(D—kV o.  Therefore,  ail  three  resonances  are  satisfied 
simultaneously,  which  then  give  rise  to  the  resonance  be¬ 
havior  of  the  emission  power  versus  magnetic  field.  Al¬ 
though  the  background  plasma  density  was  as  high  as 
10 12  cm-3,  they  believed  that  plasma  only  provided  the 
space-charge  neutralization.  However,  experimental  pa¬ 
rameters  give  fcoFo“3.54ft*  while  the  observed  mi¬ 
crowave  frequency  is  very  close  to  the  plasma  cyclotron 
frequency.  To  satisfy  their  resonance  condition,  coupling 
to  higher  harmonics  of  ko  and  flp  are  required,  which  is 
unlikely.  On  the  other  hand,  Fig.  1(a)  indicates  that 
with  sufficient  initial  transverse  momentum  the  fast- 
beam  cyclotron  mode  should  be  able  to  excite  the  fast- 
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FIG.  4.  Fast-beam  cyclotron  wave  instability  (fl*  “0.66<uf ): 

(a)  The  time  evolution  of  wave  energy  of  most  unstable  modes; 

(b)  the  time  evolution  of  beam  kinetic  energy. 

wave  branch  of  the  TM  waveguide  mode.  To  illustrate 
this  interaction,  a  simulation  run  using  px  ”0.3, 
ft»”0.66<uf,  and  without  background  plasma  was  car¬ 
ried  out.  The  wave  energy  time  evolutions  of  the  most 
unstable  modes  are  shown  in  Fig.  4(a).  The  excited 
waves  exponentiate  at  the  growth  rate  of  to,  ”0.029cof, 
which  is  weaker  than  the  plasma  cyclotron  wave  instabil¬ 
ity  and  oscillate  at  to, ■■  l.3<of,  which  satisfies  the  fast- 
beam  cyclotron  wave-dispersion  relation  (to  m,kVo+  (lb). 
The  saturation  efficiency  attained  about  3%.  Figure  4(b) 
plots  the  time  evolutions  of  beam  transverse  and  longitu¬ 
dinal  energies  which  demonstrate  that  the  instability 
only  depletes  the  transverse  component  of  beam  energy. 
This  energy  depletion  characteristic  further  supports  that 
the  interaction  is  not  due  to  the  vacuum  BWO  mecha¬ 
nism. 

As  is  evident  from  computer  simulations,  the  excita¬ 
tion  of  plasma  cyclotron  waves  requires  no  initial  beam 
transverse  momentum.  To  reconcile  the  experimental 
observation  that  the  emission  was  peaked  at  /3j.”0.3,  a 
series  of  simulations  varying  px  has  been  carried  out. 
The  change  in  px  results  in  a  change  in  the  slope  of  the 


slow  space-charge  wave  resonance  line  (Vq)  and  conse¬ 
quently  affects  the  resonance  magnetic  field.  The  simu¬ 
lation  results  show  that  the  emission  is  peaked  at  px 
“0.2.  The  increase  in  efficiency  from  px  “0  to  px  ”0.2 
is  relatively  mild  but  the  decrease  from  px  ”0.2  to 
/Jj.™ 0.3  is  more  drastic.  The  emission  is  now  peaked  at 
Clp  ”  1 ,26<uc  for  px  ”0.2  and  the  fast-beam  cyclotron  in¬ 
stability  is  overwhelmed  by  the  plasma  cyclotron  insta¬ 
bility.  In  any  event,  the  accuracy  that  the  px  measure¬ 
ment  can  attain  is  unknown  and  so  the  dependence  of  the 
emission  on  px  will  not  be  pursued  further. 

In  conclusion,  through  computer  simulations  plasma 
cyclotron  waves  are  observed  to  be  readily  excited  by  the 
relativistic  electron  beam  propagating  in  a  plasma-filled 
BWO.  When  the  excited  mode  also  coincides  with  the 
electromagnetic-wave-dispersion  relation  of  a  rippled- 
wall  waveguide,  the  slow  plasma  cyclotron  wave  is  effi¬ 
ciently  converted  into  a  fast  wave  through  the  rippled- 
wall  boundary  and  can  be  easily  detected  outside  the 
plasma.  This  mechanism  can  be  utilized  to  explain  the 
early  experimental  results  which  displayed  strong  depen¬ 
dence  of  the  emitted  power  on  magnetic  field  from  a 
plasma-filled  BWO.  Because  of  its  high  gain,  the  ob¬ 
served  mechanism  can  be  employed  to  build  compact  mi¬ 
crowave  tubes  with  medium  beam  energy  of  the  order  of 
100  kW.  The  limitation  that  the  observed  instability  will 
impose  on  the  plasma-filled  gyrotron  should  also  be  in¬ 
vestigated. 
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Abstract 


Computer  simulations  of  a  100  GHz  electron  cyclotron  autoresonance  maser 
(CARM)  oscillator  that  has  been  designed  for  an  experiment  at  the  Naval 
Research  Laboratory  are  carried  out  using  a  time -dependent,  3-D,  particle-in¬ 
cell  code.  The  results  indicate  an  intrinsically  low  efficiency  for  the 
experimental  parameters.  Enhancement  of  the  efficiency  by  varying  the  beam 
current,  interaction  length,  magnetic  field,  and  tapering  the  magnetic  field 
are  investigated.  It  is  demonstrated  that  an  efficiency  above  20%  may  be 
obtained  by  reducing  or  tapering  the  magnetic  field. 
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Theoretical  and  numerical  studies  of  the  electron  cyclotron  autoresonance 
maser  (CARM)  indicate  that  it  poses  a  strong  alternative  to  the  free  electron 
laser  as  a  high  power,  high  efficiency  radiation  source  in  the  millimeter  to 
submillimeter  wavelengths  [1-5].  However,  there  have  been  few  CARM 
experiments  to  date  and  the  potential  high  efficiency  of  the  CARM  has  yet  to 
be  demonstrated  in  the  laboratory.  The  first  CARM  experiments  [6,7]  achieved 
only  2-4%  efficiency.  More  recently,  attempts  have  been  made  at  the  Naval 
Research  Laboratory  [8]  to  develop  a  100  GHz,  24  MU  CARM  oscillator  operating 
with  20%  efficiency.  Initial  results  show  disappointingly  low  efficiency  of 
only  a  few  percent,  with  mode  competition  suspected  to  be  the  cause  [9]. 

We  have  developed  a  3-D,  time -dependent,  particle- in-cell  code  to  model  a 
CARM  oscillator  with  the  same  configuration  as  the  NRL  experiment  and  carried 
out  simulations  for  similar  beam  parameters.  The  parameters  of  the  NRL 
experiment  are  shown  in  Table  1.  The  CARM  oscillator  consists  of  a 
cylindrical  waveguide  with  a  smooth-wall  central  section  connected  on  either 
side  to  rippled-wall  sections,  or  Bragg  reflectors.  The  ripple  period  is 
half  the  wavelength  of  the  desired  cavity  mode  so  that  the  latter  is  strongly 
reflected.  The  virtue  of  this  Bragg  cavity  is  that  it  confines  the  radiation 
while  allowing  unimpeded  passage  to  the  electron  beam. 

For  the  parameters  in  Table  1,  the  cyclotron  beam  line  intersects  the 
vacuum  TEg^  mode  at  two  frequencies,  as  shown  in  Fig.  1;  the  higher  frequency 
is  the  desired  CARM  operation,  and  the  lower  frequency  corresponds  to  the 
gyrotron  interaction.  The  simulations  here  include  multiple  longitudinal 
modes  but  only  one  transverse  waveguide  mode,  the  anticipated  TEg^  mode. 
Thus,  mode  competition  between  the  TE,,  gyrotron  and  CARM  modes  is  included 
self-consistently  but  mode  competition  between  different  transverse  modes  is 


not  considered. 
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In  this  paper,  we  describe  the  numerical  model  and  present  results  from 
the  simulations.  The  main  findings  are:  Even  in  the  absence  of  mode 
competition  the  efficiency  at  the  experimental  parameters  is  very  low,  only 
about  3%.  Higher  efficiency  may  be  achieved  by  either  using  a  lower  magnetic 
field  or  tapering  the  magnetic  field.  By  tapering  the  magnetic  field, 
simulations  using  a  higher  beam  current  of  400  amps  show  it  is  possible  to 
achieve  an  efficiency  above  20%  for  an  electron  beam  with  no  velocity  spread, 
and  above  10%  with  the  anticipated  maximum  velocity  spread  of  3%. 

We  briefly  describe  here  the  numerical  algorithm  used  in  modeling  the 
CARM  oscillator.  This  numerical  algorithm  has  used  successfully  in  the  past 
to  model  gyrotrons,  and  further  details  may  be  found  in  Refs.  [10-12].  In 
the  open  waveguide  structure  of  the  Bragg  cavity,  the  transverse  dependence 
of  the  fields  at  every  position  z  along  the  axis  can  be  expanded  into  a 
complete  set  of  the  local  orthogonal  transverse  electric  and  magnetic  modes 
of  the  empty  waveguide.  For  simplicity,  however,  we  assume  that  the  desired 
TE^  mode  is  the  dominant  mode  and  the  fields  are  well  represented  by 
retaining  only  that  mode  in  the  expansion.  For  a  given  TE^  mode,  the 
electric  and  magnetic  fields  may  be  written  as 
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Here,  and  a(z)  are  the  mean  and  local  wall  radii,  respectively,  and  J '  is 

the  rath  order  Bessel  function,  and  k  the  nth  root  of  J' (k  )  -  0.  The  wall 

an  a  an 

radius  can  be  an  arbitrary  continuous  function  of  axial  position,  allowing 
the  nonuniform  cavity  wall  of  the  oscillator  to  be  modeled. 

Neglecting  the  effects  of  mode  conversion  through  changes  in  the  local 
cavity  radius,  the  time  evolution  of  the  rf  fields  from  Maxwell's  equations 
is  given  by 
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Iv  being  the  total  beam  current  in  amperes. 

The  electron  beam  is  modeled  by  thousands  of  macroelectrons  whose 
trajectories  are  followed  in  time  and  space  through  the  relativistic 
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equations  of  motion  under  the  influence  of  the  rf  and  externally  applied 
magnetic  fields: 


I<SJ>  *  h 


x  B(r  )/c  , 


(4a) 


dt 


(4b) 


The  updated  macroelectron  velocity  and  position  are  then  used  to  evaluate 

the  source  current  for  the  the  rf  fields  (Eq.3a).  In  order  to  avoid  the 

singularity  at  the  cavity  axis  in  advancing  the  electron  trajectories,  the 

calculations  are  carried  out  using  Cartesian  coordinates.  The  x  and  y 

components  of  the  transverse  field  functions  e*3  and  e ^  are  given  by 
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The  longitudinal  dependence  of  the  rf  field  is  represented  by  a  Fourier 
series  expansion.  Outside  the  cavity  region,  the  em  waves  are  numerically 
absorbed.  The  polarization  of  the  modes  is  determined  self-consistently  in 
the  simulation  through  interaction  with  the  electron  beam 

The  ripples  of  the  Bragg  reflectors  are  assumed  to  be  sinusoidal  with  the 
periodicity  and  amplitudes  given  in  Table  1.  The  reflectivities  of  the  Bragg 
reflectors  determined  numerically  by  placing  a  100  GHz,  TEg^  antenna  a 
distance  away  from  each  reflector  in  an  empty  waveguide  are  99%  for  the 
upstream  and  89%  for  the  downstream  reflectors,  in  good  agreement  with  the 


6 


theoretical  predictions. 

In  the  simulations,  the  beam  electrons  are  continuously  injected  into  the 
input  end  of  the  cavity  with  prescribed  parameters,  and  removed  at  the  output 
of  the  cavity;  the  input  beam  voltage  is  constant.  An  annular  electron  beam 
is  assumed,  with  guiding  center  located  at  0.8rw  in  order  to  maximize  the 
coupling  between  electrons  and  the  TEg^  waveguide  mode. 

The  initial  rf  field  in  the  cavity  arises  from  numerical  noise  in  the 
electron  beam.  However,  in  order  to  limit  the  computational  time  required 
to  follow  the  oscillation  to  a  steady  state  to  a  managable  value,  we  also 
prime  the  oscillator  by  injecting  an  rf  signal  into  the  cavity  for  a  short 
period  of  time  along  with  the  electron  beam.  In  all  the  simulations 
presented  below,  a  100  GHz,  -  200  kW  wave  was  injected  into  the  cavity  for 
-0.8  to  -1.8  ns  to  initiate  the  CARM  interaction.  Also,  except  for  the 
results  in  Table  6,  the  electron  beam  has  no  velocity  spread  in  the 
simulations;  all  other  parameters  are  as  given  in  Table  1,  unless  where 
otherwise  specified. 

Figure  2a  illustrates  a  cross-sectional  view  of  the  electron  beam  in  the 

cavity,  and  2b,  electron  energies  along  the  axis  of  the  cavity.  The 

electrons  are  injected  at  the  entrance  to  the  Bragg  cavity,  z  -  0,  and 

removed  at  the  exit,  z  -  7.1  cm.  Figure  3a  shows  an  example  of  the  electric 

field  profile  in  the  cavity,  3b  the  time  history  of  the  rf  field  energy,  and 

2c,  energy  of  the  spent  electron  beam  as  a  function  of  time.  The  time  taken 

to  reach  saturation  is  about  25  ns.  The  initial  sharp  rise  in  rf  field 

energy  in  Fig.  3b  (and  the  initial  drop  in  Fig.  3c)  is  due  to  the  injected 

signal.  The  efficiency  reached  is  only  about  3%.  The  efficiency,  rj ,  is 

defined  as  ij  -  (7  -7  )/(y  -1)  where  7  and  7  are  the  initial  and  spent 

o  s  o  os 

electron  beam  Lorentz  factors. 
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We  carried  out  a  parameter  search  for  better  efficiency.  Table  2  shows 
the  efficiency  and  the  (exponential)  growth  rate  of  the  rf  field  energy  in 
the  cavity  as  a  function  of  beam  current,  and  Table  3,  as  a  function  of 
cavity  length  (changing  only  the  length  of  the  smooth  section) .  None  of 
these  variation  has  much  impact  on  the  efficiency,  although  increasing  the 
current,  or  cavity  length,  increases  the  growth  rate  as  would  be  expected. 

However,  the  efficiency  could  be  increased  significantly  by  tapering,  or 
detuning  (reducing)  the  magnetic  field.  Since  tapering  or  detuning  generally 
increases  the  time  taken  to  reach  a  steady  state,  the  following  studies  were 
carried  out  using  an  electron  beam  current  of  400  A  instead  of  200  A,  and 
thus,  a  higher  growrh  rate,  in  order  to  make  the  simulation  more  tractable. 
Table  4  shows  the  efficiency  when  the  external  magnetic  field  is  uniformly 
reduced,  and  Table  5  shows  the  efficiency  as  a  function  of  the  tapered 
magnetic  field.  The  magnetic  field  is  tapered  by  reducing  it  linearly  from  z 
-  0.0  to  z  -  7.1  cm;  in  Table  5,  AB/Bq*  (B(z-7.1cm)  -  B(z-0.0))/B(z-0.0)  .  It 
is  seen  that  it  is  possible  to  achieve  efficiencies  somewhat  above  20%.  The 
sharp  drop  in  efficiency  in  Table  4,  from  21.5%  to  0%  over  1  kG  is  consistent 
with  previous  simulation  results  of  magnetic  field  detuning  [13];  although 
the  efficiency  increases  as  the  magnetic  field  is  decreased,  the  growth  rate 
decreases  and  eventually  vanishes.  Note  that  changing  the  magnetic  field 
would  change  the  frequency  in  a  superradiant  CARM  amplifier  but  in  this 
oscillator,  the  operating  frequency  is  fixed  by  the  periodicity  of  the  Bragg 
reflectors  and  is  not  affected  by  the  magnetic  field;  thus,  one  still  obtains 
100  GHz  radiation  at  the  higher  efficieny  although  the  cyclotron  resonance 
beam  line  has  shifted. 

Table  5  shows  the  effect  of  having  axial  velocity  spreads  of  1%  and  3% 
(but  with  no  energy  spread)  on  the  best  case  of  23%  efficiency  obtained  by 
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Capering  the  magnetic  field.  A  1%  spread  reduces  the  growth  rate  somewhat 
but  has  no  effect  on  the  final  efficiency.  With  the  anticipated  maximum 
velocity  spread  of  3%,  the  efficiency  drops  to  11%.  Thus,  with  a  400  A 
electron  beam,  one  might  still  be  able  to  achieve  above  20MW  of  power  within 
the  anticipated  range  of  velocity  spread. 

No  gyrotron  mode  competition  was  observed  in  any  of  the  above 
simulations.  The  efficiencies  were  thus,  obtained  under  essentially  ideal 
single-mode  operation.  The  absence  of  gyrotron  mode  competition  could  be  due 
to  our  priming  the  oscillator  at  the  CARM  frequency,  which  allows  the  CARM 
mode  to  start  at  a  much  higher  power  relative  to  the  gyrotron  mode  (which  has 
to  start  from  the  initial  noise).  We  find  that  without  the  priming,  the 
initial  broad  band  kz*  spectrum  of  the  rf  field  from  noise  will  evolve  into 
two  peaks  centered  at  the  CARM  and  gyrotron  resonances,  and  the  rf  field 
takes  a  long  time  to  grow.  When  the  simulation  is  repeated  with  an 
artificial  k  -dependent  resisitivity  to  suppress  the  gyrotron  mode,  the 
growth  of  the  rf  field  is  found  to  be  enhanced  which  suggests  that  the  the 
delay  in  oscillation  is  due  to  mode  competition.  However,  one  should  be 
cautious  about  extrapolating  this  observation  to  the  laboratory:  numerical 
simulations  typically  produce  a  more  noisy  environment  than  is  actually  the 
case  in  the  laboratory  (thus,  allowing  the  gyrotron  mode  to  start  with  higher 
power) ,  and  so  tends  to  exaggerate  the  potential  for  mode  competition. 
However,  if  priming  is  the  reason  for  the  absence  of  the  gyrotron  mode,  it 
could  be  easily  included  in  the  laboratory.  We  shall  investigate  this  issue 
in  the  future. 

In  summary,  we  have  carried  out  time -dependent,  PIC  computer  simulations 
of  a  CARM  in  a  realistic  3-D  oscillator  configuration.  Simulations  of  the 
NRL  100  GHz  CARM  oscillator  experiment  suggest  that  the  efficiency  at  the 
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experimental  parameters  is  intrinsically  low;  even  when  no  mode  competiton  is 
present,  the  efficiency  is  only  about  3%.  The  most  effective  way  to  achieve 
higher  efficiency  appears  to  be  to  reduce  or  taper  the  magnetic  field.  In 
particular,  it  is  shown  the  oscillator  could  produce  24  MW  power  with  a  3% 
beam  velocity  spread  by  using  a  beam  current  of  400  A  and  tapering  the 
magnetic  field  down  by  20%  along  the  length  of  the  cavity.  The  absence  of 
gyrotron  mode  competition  when  the  oscillator  is  primed  with  the  desired  CARM 
mode  raises  the  possiblity  that  priming  could  be  beneficial  in  the  laboratory 
in  discriminating  against  unwanted  longitudinal  as  well  as  transverse  modes 
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Table  1:  Experimental  Parameters 


Beam  voltage  [kV] 

Beam  current  [A] 

Axial  velocity  spread 

Magnetic  field  B  [kG] 
o 

Operating  mode 


600 

200 


Av  /v  <  3% 

2  Z 


25 

TE 

0.6 


61 


Cavity  parameters: 

Mean  wall  diameter  [cm]  1.59 

Upstream  reflector 

Length  [cm]  3.0 

Ripple  depth  [mm]  0.25 

Ripple  period  [mm]  1.68 

Reflectivity  99% 

Downstream  reflector 

Length  [cm]  1.5 

Ripple  depth  [mm]  0.31 

Ripple  period  [mm]  1.68 

Reflectivity  90% 

Center  Section  Length  [cm]  2.6 
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Table  2:  Efficiency  dependence  on  current 


Current  (A) 

200 

300 

600 

1000 

2000 

Efficiency  [%] 

3.3 

5.5 

6.0 

6.0 

5.0 

Growth  rate  [MHz] 

40 

184 

313 

558 

757 

Table  3:  Efficiency  dependence 

on  interaction  length 

Length  [cm] 

7.1 

9.7 

12.3 

Efficiency  [%] 

3.3 

6.0 

5.0 

Growth  rate  [MHz] 

40 

290 

422 
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Table  4:  Effieciency  dependence  on  magnetic  field  detuning 


Change  In  B-field  [kG]  0.0 

-2.0  -3.0  -4.0 

Efficiency  [%]  5.5 

15.5  21.5  0.0 

Table  5:  Efficiency  dependence  on 

magnetic  field  tapering 

AB/B  0 . 0 

o 

-0.1  -0.2  -0.25 

Efficiency  [%]  5.5 

18.0  23.0  14.0 

Table  6:  Efficiency  dependence  on  velocity  spread 

velocity  spread  [%]  0  1 

3 

Efficiency  [%]  23.0 


23.0 


11.0 
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Figure  Captions 

Figure  1.  Dispersion  relation  for  the  parameters  in  Table  1,  showing  the 

TEg^  and  cyclotron  beam  lines. 

Figure  2.  (a)  Cross  sectional  view  of  electron  beam;  is  the  cavity 

radius . 

(b)  Electron  energies  along  the  axis  of  cavity  at  steady  state 
operation. 

Figure  3.  (a)  Example  of  electric  field  profile  along  the  cavity 

(arbitrary  units) . 

(b)  Temporal  evolution  of  rf  field  energy,  U(t)  (arbitrary 
units),  in  cavity. 

(c)  Temporal  dependence  of  spent  electron  beam  energy. 


z  (cm) 


FIGURE  2 
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Abstract 

A  unified  small-signal  amplification  theory  is  developed  to  compare 
growth  mechanisms  responsible  for  a  number  of  relativistic  radiation  gen¬ 
erators.  The  theory  is  formulated  from  the  basis  that  the  electron  resonance 
frequency  produced  by  the  external  fields  of  the  devices  depends  on  ■y~q, 
where  7  is  the  beam  Lozentz  factor  and  g  is  a  constant  (g  =  1  for  cy¬ 
clotron  autoresonance  masers,  g  =  1/2  for  ion-channel  lasers  and  g  =  0 
for  free  electron  lasers).  It  is  concluded  that  for  wave  amplification,  the 
sign  of  the  electron  mismatch  frequency  is  required  to  be  the  same  as  the 
sign  of  bunching  parameter  which  is  determined  by  the  total  bunching 
both  axial  and  azimuthal;  this  depends  on  the  q  value.  The  two  bunching 
mechanisms  exist,  not  only  in  the  single  electron  resonance  regime,  but 
also  in  the  collective  gain  regime.  Competition  or  reinforcement  between 
the  two  bunching  mechanisms  is  determined  by  the  q  value,  the  electron 
axial  velocity  and  the  wave  phase  velocity. 


PACS  numbers:  42.52,+x,  42.55.Tb,  52.40.Mj,  52.75.Ms 


1.  Introduction 


Electron  bunching  is  the  amplification  mechanism  responsible  for  many  state  of  the 
art  devices  for  generating  radiation  from  microwaves  to  the  ultraviolet.  Among  the 
many  available  devices,  cyclotron  masers  (CARMs  and  gyrotrons)  and  free  electron 
lasers  (FELs)  are  the  most  successful. 

Electron  cyclotron  maser  instabilities,  independently  proposed  by  Twiss[l], 
Schneider[2]  and  Gaponov[3]  in  the  late  50’s,  involve  electrons  gyrating  about  an 
external  axial  magnetic  field.  At  the  same  time  Weibel[4]  studied  electron  whistler 
instabilities  driven  by  an  anisotropic  velocity  distribution  of  electrons.  Both  involve 
the  transverse  electromagnetic  wave  fields  interacting  with  gyrating  electrons.  Two 
decades  later,  Chu  and  Hirshfield(5]  presented  a  unified  dispersion  relation  and  pointed 
out  that  net  bunching  resulted  from  a  competition  between  azimuthal  (bunching  due  to 
cyclotron  frequencies  changes  associated  with  the  relativistic  mass  effect)  bunching  and 
axial  (conventional  overtaking)  bunching  is  responsible  for  these  instabilities.  Bratman 
et  al.[6]  and  Lin[8]  investigated  the  electron  cyclotron  maser  with  emphasis  on  the 
Doppler  shift  effect  (the  so  called  cyclotron  autoresonance  maser  or  CARM). 

In  the  late  50’s  and  early  60’s,  Motz  and  Nakamura[13]  and  Phillips[14]  studied 
the  radiation  mechanism  for  an  electron  beam  going  through  a  periodic  undulators.  In 
1971  Madey[15]  described  the  stimulated  emission  mechanism  of  a  free  electron  laser 
based  on  a  relativistic  electron  beam  moving  through  a  periodic  transverse  dc  magnetic 
field  and  later  carried  out  successful  experiments  in  1976[16]  and  1977[17]  at  Stanford. 
Colson[18]  presented  a  one-body  analysis  of  the  FEL  by  solving  the  pendulum  equation. 
Dispersion  relations  have  been  derived  for  the  collective  regime  of  the  FEL  by  Sprangle 
et  al.[19].  Computer  simulations  by  Kwan  et  al.[20],  Lin  and  Dawson[21]  and  Liewer 
et  al.[22]  have  been  used  to  check  linear  theory,  nonlinear  saturation  mechanisms  and 
effeciencies  of  FELs. 
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CARMs[l]-[12]  and  FELs[13]-[24]  have  independendy  been  the  subjects  of  con¬ 
tinuous  activity  for  the  past  three  decades.  However,  the  similarities  and  differences 
between  the  two  amplification  mechanisms  involved  remain  to  be  further  clarified. 

Recently  the  Ion-Channel  Laser  (ICL)  as  another  means  to  generate  tunable  coherent 
radiation  was  conceived[25,26];  in  this  device  a  relativistic  electron  beam  is  injected 
into  an  underdense  plasma  (the  plasma  density  is  less  than  or  equal  to  the  beam  density). 
The  beam  front  pushes  out  the  plasma  electrons  leaving  an  ion  channel.  Ion  focusing 
causes  the  beam  electrons  to  oscillate  about  the  axis  and  plays  a  similar  role  as  the 
magnetic  field  in  a  CARM.  Using  this  type  of  device  it  appears  that  radiation  can  be 
produced  with  wave  frequencies  ranging  from  microwaves  to  soft  X-rays;  different 
devices  must  be  used  for  different  frequency  ranges. 

This  study  is  motivated  by  growing  activity  on  relativistic  devices  and  was  also 
triggered  by  findings  on  the  amplification  mechanism  of  ICLs.  The  goal  is  to  develop 
a  unified  small-signal  gain  theory  to  illuminate  the  physical  issues  associated  with  the 
amplification  mechanisms  of  CARMs,  FELs  and  ICLs,  and  to  place  them  on  a  common 
ground.  The  electron  beams  being  studied  are  mono-energetic,  have  no  velocity  spread 
and  gyrate  along  their  drifting  direction.  The  electromagnetic  wave  is  assumed  to  be 
right  circularly  polarized  and  its  amplitude  is  small  enough  to  validate  perturbation 
theory.  The  rate  of  electron  energy  change  and  the  wave  energy  gain  are  evaluated  to 
illustrate  the  amplification  mechanism. 

CARMs,  ICLs  and  FELs  all  have  a  zeroth  order  radial  force  on  the  electrons  which 
induces  oscillatory  motion  about  their  overall  drift.  The  resulting  electron  frequencies 
have  different  dependencies  on  the  electron  relativistic  factor,  depending  on  7"’.  The 
value  of  exponents  (-g)  of  7  reveals  different  degrees  of  azimuthal  bunchings  caused 
by  the  electron  relativistic  mass  variation.  Also  FELs  differ  from  CARMs  and  ICLs  in 
that  they  have  an  additional  axial  pondermotive  force  induced  by  the  magnetic  wiggler. 

Some  interesting  physics  insights  will  be  discussed.  Both  inertial  bunching  and 
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force  bunching[6,7]  are  included  in  our  calculation.  Force  bunching  can  be  neglected 
under  some  conditions.  The  inertial  bunching  has  two  parts  :  azimuthal  bunching  and 
axial  bunching.  The  axial  force  will  not  give  rise  to  electron  axial  density  bunching  if 
initially  the  electrons  have  uniform  azimuthal  phases.  As  a  result  there  will  be  no  space 
charge  force.  The  azimuthal  bunching  which  relies  on  relativistic  effect  is  induced 
by  the  wave  electric  field.  On  the  other  hand,  the  axial  bunching  is  a  result  of  the 
competition  between  the  conventional  v  x  B  force  and  the  axial  velocity  change  caused 
by  the  relativistic  mass  variation.  For  CARMs  and  ICLs,  the  relativistic  mass  variation 
effect  on  the  axial  velocity  is  dominant  if 

c2  <  v0zvph  (1) 

where  c  is  the  speed  of  light,  t»oz  is  the  beam  axial  velocity  and  vph  is  the  phase  velocity 
of  wave.  This  is  because  the  beam  velocity  is  greater  than  the  wave  group  velocity. 
We  will  show  that  for  FELs  the  conventional  v  x  B  force  is  always  dominant  due  to 
the  pondermotive  force  of  a  magnetic  wiggler. 

The  competition  of  different  bunching  mechanisms  can  be  seen  not  only  in  the 
single  electron  resonance  condition^],  but  also  in  the  rate  of  electron  energy  change. 
Azimuthal  bunching  and  axial  bunching  do  not  always  offset  each  other  as  concluded 
in[5].  If  condition  (1)  is  satisfied,  they  may  reinforce  each  other.  The  competition  will 
be  shown  to  depend  on  the  beam  axial  velocity  vo,  (and/or  energy)  for  ICLs,  but  is 
independent  of  vo*  for  CARMs. 

By  evaluating  the  electron  phase  change[9,27]  and  using  perturbation  theory,  we 
have  developed  a  unified  small-signal  theory  which  gives  a  same  formulation  for 
CARMs,  ICLs  and  FELs  and  reveals  the  general  characteristics  for  relativistic  radi¬ 
ation  generation.  The  first  order  results  contain  the  physical  meaning  of  azimuthal  and 
axial  bunchings  while  the  rate  of  electron  energy  change  and  the  wave  energy  gain  can 
be  obtained  from  second  order  results  (<  v\  E\  >)  which  also  show  the  competition  of 
the  different  bunching  mechanisms.  The  same  result  can  also  be  obtained  by  evaluating 
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the  phase  shift  of  an  electron  relative  to  its  resonance  phase  with  a  wave.  In  Section  III., 
the  similarities  and  differences  among  CARMs,  ICLs  and  FELs  are  illustrated  through 
our  unified  theory.  Section  IV.  gives  a  summary  and  discussion. 
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11.  Unified  Small-Signal  Theory  for  Relativistic  Radiation 


Consider  an  electron  undergoing  transverse  oscillation  as  it  moves  along  the  axial 
direction.  There  is  a  zeroth  order  force  in  the  radial  direction  which  may  come  from 
an  axial  magnetic  field  Bo  as  in  CARMs,  a  transverse  electric  field  Ex  as  in  ICLs  or  a 

wiggler  magnetic  field  Bw  as  in  FELs.  For  simplicity,  assume  the  force  is  centripetal 

which  results  in  the  electron  motion  having  a  helical  trajectory.  As  shown  in  Fig.l.  the 
radial  force  and  electron  velocity  can  be  expressed  by 

F  =  F[-£  cos  <f>  —  y  sin  4>]  (2) 

v  =  fj.[— £sin<£  + j/cos<£]  +  u*2  (3) 

where  4>  =  +  4>0  is  the  electron  oscillatory  (azimuthal)  phase  angle,  =  4>{t  =  0), 

u,  -  U'/y  is  the  electron  oscillating  frequency  which  is  dependent  on  its  Lorentz  factor 
7  =  (1  —  v^/c2— v2z/c2)~l  with  exponent  -q,  is  the  non-relativistic  electron  oscillating 
frequency,  F  =  is  the  amplitude  of  the  radial  force,  me  is  the  electron  rest 

mass,  vx  is  the  amplitude  of  the  electron  transverse  velocity,  and  vz  is  the  electron  axial 
velocity.  Following  common  assumption  we  assume  that  the  electron  beam  is  tenuous, 
thus  electron  space  charge  is  not  important  and  the  role  of  electrostatic  waves[10]  can 
be  ignored  in  the  amplification  mechanism.  We  consider  a  right  circularly  polarized 
plane  electromagnetic  wave;  its  Doppler  shifted  frequency  is  comparable  to  the  electron 
oscillation  frequency; 

Ei  -  £i[zcos(u><  —  kz)  +  y sin(u><  —  kz)]  (4) 

B\  =  B|[— £sin(u;t  —  kz) +  y  cos  (ut  —  kz)}  (5) 

where  E\{B\)  is  the  amplitude  of  wave  electric  (magnetic)  field,  u>  is  the  wave  frequency, 
k  *  kz  is  the  wavenumber  and  z  is  the  axial  spatial  position,  ty,  =  w/fc  is  the  wave 
phase  velocity. 
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At  this  point,  we  have  not  introduced  any  ordering.  Under  these  assumptions,  the 
governing  equations  for  an  electron  interacting  with  the  wave  are 


d  , 

—(7  mec2)  = 

—eE  |Vx  sin(ui<  —  kz  —  <f>) 

(6) 

1 

rift. 
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-e£,(  1  -  —)  +  F 

(7) 

at 

±P  = 
dt  1 

U ) 

--{vl  X  B\)  +  Fp 
c 

(8) 

d 

dt 

Vt 

(9) 

where  — e  is  the  electron  charge,  P±  =  jmev~l  is  the  electron  transverse  momentum, 
Fp  =  0  for  CARMs  and  ICLs,  Fp  =  —  *(v±  x  Bw)  is  the  pondermotive  force  of  magnetic 
wiggler  for  FELs,  and  Pz  =  7 mtvz  is  the  electron  axial  momentum.  Here  we  should  be 
aware  that  the  exchange  of  energy  between  the  wave  and  electrons  is  associated  with 
changes  in  electron  velocities. 

By  using  Faraday’s  law  B\  =  ^{k  x  E\),  the  pondermotive  force  due  to  the  wave 
magnetic  field  can  be  expressed  by 

-  -(vl  x  B\)  =  -77 yni'C 2  .  (10) 

c  u>  at 


In  an  FEL,  as  v±/c  <x  7“',  and  from  the  definition  of  7  we  obtain 

%  =  .  (ii) 

at  vz  7  at 

Comparing  Eq.(ll)  with  Eqs.  (8)  and  (10)  the  pondermotive  force  of  a  magnetic  wiggler 
can  be  written  by 

Fp  =  (12) 

u  at 

where  kw  is  the  wavenumber  of  the  magnetic  wiggler. 

Combining  Eqs.(10-12)  and  Eq.(8),  the  axial  momentum  equation  becomes 


d  _  1  d  2 

-yPr  =  — — 

at  vp  at 


(13) 
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where  vp  =  u>/k  for  CARMs  and  ICLs  and  vp  =  u/(k  +  kw)  for  FELs. 

Now  we  are  going  to  do  ordering.  Since  E\  is  a  small  parameter,  e£iux/(7mec2)  << 
1,  the  governing  equations  can  be  solved  by  iteration.  To  zeroth  order,  Eqs.  (6)-(9)  and 
(13)  can  be  reduced  to 


7  =  70  (14) 

t<i  =  voj.[-zsin</>  +  ycos  <j>]  (15) 

vt  -  v0z  (16) 

<j>  =  utrft  +  fo  (17) 

z  =  V0zt  +  Zo  (18) 


where  70  =  (1  -  t>oi/c2  ~  ^Oi/c2)-*  is  the  unperturbed  electron  7,  vox(i’oz)  is  the 
amplitude  of  the  unperturbed  electron  transverse  (axial)  velocity,  u ;f0  =  we/7o  is  the 
unperturbed  electron  oscillatory  frequency  and  zq  is  the  initial  electron  axial  spatial 
position  while  is  the  initial  electron  azimuthal  phase  angle. 

To  obtain  the  first  order  solution  for  7,  we  substitute  Eqs.  (17)  and  (18)  into  Eq.  (6) 
and  obtain 

,  d 

m'C  di 7=  -e-Eivo!  sm[(w  -  kvoz  -  w#)t  -  (fcz°  +  <£o)]  .  (19) 

Substituting  Eqs.  (17),  (18)  and  (19)  into  Eq.  (13)  yields  the  first  order  solution  for 

Vt- 

To mt^rvz  =  -eE|Vox[—  -  ^]sin[(u>  -  kv 0l  -  u >/>)<  -  ( kzc  +  <£o)]  .  (20) 

dt  vp  <r 

Substituting  Eq.  (19)  into  Eq.  (7)  and  separating  the  equations  for  the  changes  of 
the  perturbed  electron  transverse  velocity  and  azimuthal  phase  angle  gives 

4*±i  =  __££L(i_^£_l4)sin[(u;-fcVol-^)<-(fczo  +  ^)]  (21) 

dt  Torn,  u>  & 

~4>  =  H ’•+  eE]  (1  _  ^2£)cos  [(u>  -  ki*t  -  u>,0 )/  -  ( kzo  +  fr,))  .  (22) 

dt  7’  7omevoi 
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(23) 


Eqs.  (19)  and  (20)  are  readily  integrated  to  give  7  and  vt  to  first  order  in  Et: 


7  =7o  + 
vz  =  VOz  + 


eEivo±  cos  [(u>  —  kvoz  —  u,o)t  —  (kzo  +  fa)]  —  cos  (kzo  +  fa) 


mr 

eEivo±.c2 

- rl - v0J 

70 me-4  Up 


w  —  kv oz  —  w,o 


cos  [(u>  -  kvpz  —  o^k))/  -  (kzp  +  </>q)]  -  COS  ( kzp  +  fa) 
u>  —  kvoz  — 


(24) 


Substituting  Eqs.  (23)  and  (24)  into  Eqs.  (22)  and  (9)  respectively  and  integrating 
Eqs.  (21)-(22)  and  (9)  yields  the  expression  for  the  perturbed  amplitude  of  the  electron 
transverse  velocity,  its  azimuthal  phase  angle  and  the  axial  phase  position: 

v$±  cos  V’  -  cos  ( kzo  +  fa) 


eE  1  _  kvQ. 

7ome  u> 


v±i 

<f)  =  fa  +  u>,o  t  —  qu>A 


c2  ^  (u;  —  kvoz 

eEwox 


^to) 


{ 


70^^ 

sin  ip  +  sin  (kzp  +  <f>o)  _ 

(u>  -  kvoz  —  u>,o)2  u>  —  kvoz  —  u>,o 


f  cos  ( kzo  +  fa) 


} 


eE\ 


^  kvoz^smxl’  +  sm(kzo  +  fa) 
lorn  t  vo  ±  u>  u>  —  kv0i  —  u>,o 


c*  .eEjUoi 


z  =  zo  +  u0lt  +  [ - voJ 


{ 


7 0mec2 

sinV»  +  sin(fczo  +  ^o)  tcos(kzp  +  <ft0) 
(u>  -  kvoz  —  *a)2  w  —  kvQz  —  u>,o 


} 


(25) 


(26) 


(27) 


where  rp  =  (u;  -  fcu0z  -  w,o)t  -  (fcz 0  +  <£o).  The  perturbed  amplitude  of  the  electron  trans¬ 
verse  velocity  uxi  and  the  fourth  term  of  Eq.  (26),  fa /  (the  azimuthal  phase  variation  due 
to  forces  acting  on  the  electron  in  the  radial  direction),  are  called  force  bunching[6,7]. 
On  the  other  hand,  the  axial  phase  variation  z\  and  the  third  term  of  Eq.  (26),  <pu  (the 
azimuthal  phase  variation  due  to  the  relativistic  mass  effect  on  the  electron  oscillatory 
frequency),  are  called  inertial  bunching  because  they  are  related  to  coordinate  frame. 
4>\j  is  generally  greater  than  v±\/voi-  From  Eq.  (26)  the  force  bunching  is  negligible 
if  v2Jc2  >  >  (u>  -  kvoz  -  u>io)/u>,  that  is,  when 


where  o>,0f  %  2ttx  the  number  of  electron  oscillation  periods.  This  inequality  is  usually 
satisfied  in  regimes  of  interest. 

For  (u i  -  kvoz  -  u),o )t  1,  which  is  true  for  resonant  interaction  at  early  time,  if 

we  neglect  the  force  bunching  term,  Eq.  (26)  and  Eq.  (27)  can  be  further  simplified  to 

eE\  vox  f2 

<fr  ~  (fro  +  u>,ot  q^io - =•— sin(kzo  +  <£o)  (29) 

7o  me<r  2 

c2  e  Ei  vox  f2 

z  ~  zo  +  votf  +  [ - vo2] - sin(fczo  +  (fro)  •  (30) 

vp  7ome<r  2 

Although  the  first  order  solution  gives  no  net  energy  exchange  between  the  wave  and 
electrons  ar  can  be  seen  by  averaging  Eq.  (23)  over  the  initial  phases  zo  and  <fro  of  all  the 
electrons,  it  illustrates  that  the  electrons  are  bunched  in  azimuthal  phase.  Eq.  (29)  reveals 
that  at  any  position  (e.g.,  zo  =  0)  electrons,  initially  from  the  upper  half  plane  (0  < 
(fro  <  7r),  slip  back  in  phase  due  to  the  mass  increase  imparted  by  the  wave  eleuric  force 
while  the  electrons  initially  from  the  lower  half  plane  (?r  <  (fro  <  27r)  advance  in  phase 
due  to  the  mass  decrease;  thus  the  electrons  are  bunched  in  azimuthal  phase.  Azimuthal 
bunching  can  occur  in  electron  beams  initially  with  uniform  or  coherent  azimuthal  phase 
distribution.  As  implied  by  Eq.  (30),  only  for  nonuniform  phase  angle  plasmas  does  the 
axial  force  produce  axial  spatial  bunching.  On  the  other  hand,  the  axial  force  will  not 
give  rise  to  a  spatial  bunching  for  uniform  phase  angle  plasmas.  This  is  because  after 
taking  the  initial  azimuthal  phase  angle  average  the  density  perturbation  becomes  zero 
even  though  every  individual  electron  has  been  spatially  perturbed.  However  the  change 
in  electron  axial  position  will  result  in  azimuthal  phase  bunching  with  respect  to  the 
wave  through  the  Doppler  effect.  Since  energy  exchange  involves  phase  perturbation, 
both  azimuthal  (relativistic)  bunching  and  axial  (conventional)  bunching  can  play  an 
important  role  in  determining  the  rate  of  electron  energy  change  to  second  order. 

To  second  order  the  energy  change  is 

•j-y mtc2  -  — eE|Vox  sint/» 
at 

— eEi(— voikzt  cosifr  —  vox^i  cos xfr  +  vXi  sin  V’l  (31) 
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where  <p\  =  <pu  +  <p\j. 

Substituting  Eqs.  (25)-(27)  into  Eq.  (31),  we  obtain  the  rate  of  electron  energy 
change  to  second  order  in  E\ : 


a  2  nr 

— 7mec  =  -eEwox  { 


sin  ip 
eE\vpx 

lomeC2 


[qu, 0  -  k{ - t’o,)] 


[ 


sin  ip  cos  t/’  +  sin  ( kzo  +  <po)  cos  ip 
(u  -  kv0t  -  W/l)2 
t  COS  (fc  2o  +  <Po)  cos  ip 

(W  -  kV0z  -  Urf) 

_  kVt) 

7ome  <*> 

sin  ip  cos  ip  +  sin  {kzo  +  <ft>)  cos  V» 


e£, 


(u»  -  *v0i  -  u'.o) 
kvz 


(1  - 


l’0l 


) 


70  me  w  c2 

sin  V’  cos  V’  -  cos  {kzo  +  <£o)  sin  ip 
(w  -  fcv0l  -  u>,o) 


(32) 


Eq.  (32)  is  the  expression  for  the  energy  exchange  rate  for  a  single  electron  whose  initial 
azimuthal  (axial)  phase  angle  (position)  is  at  <po  {zq).  To  obtain  the  rate  of  net  electron 
energy  change,  we  integrate  Eq.  (32)  over  <po  and  kzo  '• 


'EM 


01 


270  mec2 


[90;/)  -  k{ - vo*)] 


[ 


sin(w  —  kv o,  —  t  cos(u>  —  fcvo*  —  tv*o)^ 


(<*>  -  kvoz  -  wa)2 


uj  —  kvoz  -  u>j) 


] 


J&m  1  -  ^  )  _  ^0XjS^n(c,;  ~  ^’0*  ~ 

27ome  u)  c2  (u>  —  kvoz  —  u?,o) 


(33) 


Wave  is  saturated  when  <  jt‘ymtc2  >=  0.  Thus,  the  saturation  time  T,at  ~  4.4/ (a^  - 
kvot  -  u>(o)  can  be  estimated  by  setting  the  inertial  bunching  term  (the  first  term  of  the 
right  hand  side  of  Eq.  (33) )  to  zero  if  we  neglect  the  force  bunching  term.  The  second 
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term  is  caused  by  force  bunching  and  vanishes  for  7ovo±/c  ~  1.  (w  -  kv o*  —  u;^)t  is 
approximately  n  at  saturation,  so  taking  (u>  -  kv oz  —  u>, o)i  of  order  1  we  find  that  the 
force  bunching  term  can  be  neglected  once  7ouo±/c2  >>  l/(w*oO-  This  is  the  same 
condition  found  previously  (Eq.  28).  From  now  on,  we  will  drop  the  force  bunching 


term.  The  time-dependent  change  of  the  mean  electron  energy  can  be  expressed  as 

e2E2v2  c2 

(7  ~  7 o)mec2  ~  - -  k( - vo*)] 

7omec2  vp 


cos  (u>  —  kvoi  —  w,o)t  —  1  +  ^ 


_  !  +  fai-fcflU.-JiMk  sin(u>  _  kv0t  -  U^)* 

7 — r — 2 — \3 - ^34> 

(u>  -  k v0t  -  u>.o Y 


From  energy  conservation,  multiplying  (70  —  7)mec2  by  the  total  number  of  electrons 


ntV  and  divided  by  the  total  wave  energy  2E2V/%tt,  where  ne  is  the  electron  density 


and  V  is  the  volume  of  interaction  regime,  yields  the  gain  expression: 

^  uLvo±r  ,  c2  cos  6t  -  1  +  Qsm  6t 

Gain  ~ - Ef~r^[quA  -  k( - vo*)J - Tv  - - 

70c2  vp  b3 


for  <  1 


*3  / 1  voz.  ttr  voz x 

B  ~  ?(1 - )  -  (-=■ - ) 


6  =  w  —  fct’Oi  — 


where  B  is  the  bunching  parameter,  b  is  the  frequency  mismatch  factor  and  uy  = 
(47rnee2/me)i  is  the  electron  plasma  frequency.  Eqs.  (35)-  (38)  explicitly  express  the 
gain  as  a  function  of  the  electron  beam  energy  70,  the  electron  azimuthal  velocity  vox. 
the  electron  plasma  frequency  w,*,  the  wave  frequency  u),  the  mismatch  6  and  the 
bunching  parameter  B.  At  the  beginning  of  the  interaction  the  gain  is  proportional  to 
the  4th  power  of  the  time.  As  the  interaction  length  is  fixed,  the  gain  is  maximum  at 
6t  ~  2.6  and  proportional  to  the  cubic  power  of  the  time. 

The  mismatch  6  =  v  —  kv 0i  -  ~  0  indicates  how  closely  the  electrons  are  to 

resonance  with  the  wave.  It  also  gives  the  wave  frequency  scaling: 
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where  yoz  =  1/(1  -  r^/c2)^  and  we  have  assumed  that  the  wave  phase  velocity  and  the 
beam  velocity  are  near  the  speed  of  light.  Eq.  (39)  shows  that  to  attain  high  frequency 
radiation  for  q  <  2  high  electron  oscillatory  frequency  or  high  energy  electron  beams 
have  to  be  employed.  But  as  indicated  in  Eq.  (35)  the  increase  in  beam  energy  results 
in  lower  gain.  Thus  a  practical  choice  for  a  high  frequency  radiation  source  is  to  use 
very  high  electron  oscillatory  frequency  and  moderate  electron  beam  energy. 

Eq.  (36)  also  indicates  that  to  achieve  wave  amplification  the  mismatch  and  the 
bunching  parameter  must  have  the  same  sign.  Therefore,  the  bunching  parameter  B 
and  mismatch  factor  6  should  be  carefully  chosen  in  order  to  give  positive  gain.  Their 
relation  is  shown  in  Table  1.  There  are  two  bunching  mechanisms  (azimuthal  and  axial 
bunchings)  involved  in  the  bunching  parameter.  The  physical  mechanism  associated 
with  the  bunching  parameter  and  mismatch  factor  can  be  also  seen  by  examining  the 
resonance  condition 

u>  —  kvz  ^  —  =  6  ~  0  .  (40) 

7* 

Instead  of  defining  the  Doppler  shift  frequency  of  wave  a s  u>  —  kvt,  for  convenience 
we  define  the  Doppler  shifted  frequency  of  an  electron  to  be 


wp  =  kvz  ±  — 


(41) 


which  is  the  electron  oscillatory  frequency  related  to  wave  in  the  beam  frame.  By 
interacting  with  the  wave  the  electron  axial  velocity  and  energy  change  which  results 
in 

*  (42) 


Aujd  =  kAvz  qp  q— ^7 
7o 

where  Au>D  =  u >£>(<  =  At)  —  u>o(t  =  0),  Avz  =  vt(t  =  At)  -  vz(t  =  0),  Ay  =  7 (t  =  At)  —  70 
and  we  have  assumed  that  the  perturbation  is  small  (i.e.,  Ay  70).  From  Eq.  (6),  the 
energy  change  is 


Ay  =  — 


mtc * 


:(t>0j.  •  E\)At 


(43) 
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From  Eq.  (6)  and  Eq.  (13),  we  obtain  the  axial  velocity  change: 

2 

^  iTil  ■  f 

Avt  =  -(-2  -  -)[-^-2voi  •  Ex) At  .  (44) 

Vp  Vp  7o  rri'C* 

Substituting  Eqs.  (43)  and  (44)  into  Eq.  (42)  and  using  the  resonance  condition  (40) 
yield 

Aud  =  B^-Av ol  •  E\ )At  (45) 

70  mecz 

which  relates  A u>d  to  the  bunching  parameter.  We  have  also  shown  (Eq.37)  that  the 
same  formula  for  the  bunching  parameter  can  be  obtained  by  using  small-signal  gain 
theory. 

The  first  term  in  the  bunching  parameter  expression  from  the  small-signal  theory 
result  (Eq.  37)  comes  from  the  variation  of  the  electron  azimuthal  phase  angle  (Eq.  26) 
induced  by  interacting  with  the  azimuthal  electric  field.  Hence,  the  first  term  of  the 
bunching  parameter  is  called  azimuthal  bunching.  On  the  other  hand,  the  second  term 
is  called  axial  bunching.  It  comes  from  the  variation  of  the  electron  axial  phase  position 
(Eq.  27)  resulting  from  the  perturbation  of  the  axial  velocity.  The  first  term  of  the  axial 
phase  variation  is  due  to  the  perturbation  of  the  axial  velocity  by  the  nonrelativistic  v  x  B 
force,  but  the  second  term  is  due  to  the  relativistic  mass  variation.  The  axial  bunching  is 
the  result  of  competition  between  the  mass  variation  effect  and  the  conventional  v  x  B 
effect. 

Figure.  2,  shows  the  relations  between  the  two  bunching  mechanisms  in  the  plot 
Vp/c  is  the  horizontal  axis  and  vqJc  as  its  vertical  axis.  It  can  be  divided  into  four 
regions.  In  region  (I),  the  beam  velocity  is  greater  than  the  phase  velocity  of  the  wave 
and  separated  from  the  others  by  the  line  vqz  =  vp.  In  this  region  azimuthal  bunching 
reinforces  axial  bunching  and  the  bunching  parameter  is  negative  (B  <  0).  In  region  (II) 
and  region  (HI)  the  two  bunchings  offset  each  other.  Axial  bunching  is  dominant  and 
the  bunching  parameter  is  negative  in  region  (II)  while  azimuthal  bunching  dominates 
and  the  bunching  parameter  is  positive  in  region  (IE).  These  two  regions  are  separated 
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by  the  curve 


/i  \Vo*  c  vp 

(1  -  q) —  = - q—  . 

c  v„  c 


As  voz/c  — *•  0,  vp/c  approaches  at  the  point  P.  On  this  curve,  the  bunching 
parameter  B  is  zero  and  there  is  no  radiation  (Eq.  36).  Using  operating  parameters 
slightly  away  from  this  curve  to  obtain  a  small  but  finite  gain,  the  electrons  can  remain 
resonant  with  the  wave  for  a  long  time  since  the  two  bunchings  effects  offset  each  other. 
This  is  called  autoresonance.  Region  (IV)  is  separated  from  region  (HI)  by  the  curve 
1  =  voz/c  •  vp/c.  On  this  curve  the  axial  bunching  term  is  zero  because  the  bunching 
from  the  v  x  B  force  is  exactly  cancelled  by  that  of  the  relativistic  mass  variation  In 
region  (IV),  the  relativistic  mass  variation  effect  on  the  axial  bunching  is  dominant 
while  the  conventional  v  x  B  force  effect  is  dominant  in  other  regions.  Hence  the 
axial  bunching  reinforces  azimuthal  bunching  and  the  bunching  parameter  is  positive 
(B  >  0). 

Although  Eqs.  (35)-(38)  and  (45)  are  derived  by  including  relativistic  mass  varia¬ 
tion,  the  derivations  are  still  valid  in  the  weakly  relativistic  limit  7  — ♦  1.  Competition 
between  the  two  bunchings  persists  not  only  in  relativistic  beams  but  also  in  weakly 
relativistic  plasmas.  The  relativistic  (azimuthal)  bunching  may  dominate  for  weakly 
relativistic  electrons  and  the  conventional  (axial)  bunching  may  dominate  at  very  rela¬ 
tivistic  electron  energies  as  implied  by  Fig.  2. 

To  obtain  the  time-dependent  mismatch  factor  we  rewrite  the  resonant  condition  (40) 


6(At)  -  u>  -  u>D(t  =  0)  -  Awd  ~  0  for  resonance 


For  wave  amplification  more  electrons  give  up  energy  ( vo±  •  E\  >0)  to  the  wave 
than  take  energy  from  it  so  that  for  most  electrons  Au >d  and  B  must  have  the  same 
sign  (Eq.  45).  In  the  case  of  positive  (negative)  B,  the  majority  of  electrons  increase 
(decrease)  their  Doppler-shifted  frequency,  that  is,  A 07?  >  0  (<  0).  Therefore,  an  initial 
positive  (negative)  mismatch  factor  S  is  essential  in  order  to  retain  most  electrons  in 
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resonance  with  the  wave  at  early  time.  This  explains  why  the  sign  of  the  mismatch  factor 
is  required  to  be  the  same  as  the  sign  of  the  bunching  parameter  for  wave  amplification 
as  shown  in  Table  I. 
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HI.  the  Comparison  of  CARMs,  ICLs  and  FELs 


In  the  previous  section,  a  unified  theory  was  developed  to  describe  the  general  wave 
amplification  process.  In  this  section,  the  unified  theory  will  be  applied  to  different 
devices.  In  a  CARM  configuration  (Fig. 3a.)  an  external  magnetic  field  ,Bo,  is  imposed 
along  the  axial  direction.  The  centipetal  radial  driving  force  which  makes  electron  gyrate 
comes  from  vox  x  Bo  force.  In  this  case  u>,  =  u>ct  where  u>ce  =  eBo/yrritC  =  QC'/l 
is  the  electron  cyclotron  frequency,  q  =  1,  and  Fq  =  eV0xBo/c  is  the  amplitude  of  the 
driving  force.  The  resonance  condition  is 

u>  -  kvz  -  —  ~  0  .  (48) 

7 

The  resonant  wave  frequency  is  a>  ~  2^zClce/y  ~  2~fQ.cr  where  we  have  assumed 
voz  >  rox-  The  rate  of  electron  energy  variation  becomes 


where  Be  is  the  CARM  bunching  parameter  and  6e  is  the  CARMs’  mismatch  factor. 
Here  the  results  are  general  and  include  the  bunching  parameter  dependence.  In  Ref.  [9], 
a  similar  result  was  obtained  except  that  the  bunching  parameter  was  equal  to  one 
which  resulted  from  assuming  vo,  =  0  and  that  the  wave  magnetic  field  B j  =  0.  Their 
calculation  missed  the  competition  between  the  two  bunching  mechanisms.  Also,  their 
coefficient  16  in  the  rate  of  energy  change  should  be  24.  This  discrepancy  arises 
because  they  neglected  the  cubic  bunching  part  (Eq.  26  -  27).  In  Ref.  [5],  the  bunching 
parameter  is  derived  by  examining  the  resonance  condition  and  the  result  was  the  same 
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as  the  second  line  of  Eq.  (50).  Their  study  concluded  that  the  two  bunchings  always 
offset  each  other  and  the  bunching  parameter  is  independent  of  the  beam  drift  velocity. 
This  can  only  apply  to  regions  (13)  and  (HI)  of  Fig.2.  In  regions  (I)  and  (TV),  the  two 
bunchings  reinforce  each  other  and  autoresonance  is  impossible.  However,  region  (II) 
and  region  (HI)  are  separated  by  the  vph  =  c  line,  which  is  independent  of  the  beam 
velocity.  Here  azimuthal  (axial)  bunching  is  dominant  in  the  fast  (slow)  wave  region. 

In  a  circular  polarized  ICL  configuration  (Fig. 3b.),  the  ion  electric  field  which  pro¬ 
vides  a  centripetal  driving  force  is  represented  by 

E,  =  A,  ■  r[cos  4>x  +  sin  <f>y]  (52) 

where  A,  =  2*nxe,  n,  is  the  ion’s  density,  r  is  the  electron  gyroradius.  A,  •  r  is 
the  amplitude  of  the  ion  electric  field  and  <$>  is  the  angular  phase  angle.  In  this  case 
Fo  =  27rnie2r  is  the  amplitude  of  the  driving  force,  v0i  =  u>,r,  q  -  {  and  -  u>ei  = 
27rnie2/me.  The  resonance  condition  is 

u>  -  kvz  —  ~  0  •  (53) 

The  resonant  wave  frequency  is  w  ~  2yhjjex.  The  rate  of  electron  energy  change  and 
the  gain  are 


e2Eh'o,v  „  rsmSit  tcos<M, 

*  2W  ®'[  «?  -  t,  ’ 

(54) 

Gain 

wptvoiw  „  2  cos  Sxt  —  2  +  6it  sin  6ti 

2jo  c2  '  S2 

(55) 

Bi 

2  Vph  Vph  Vph 

(56) 

Si 

,  ^ e » 

=  u>  —  kvoz - r 

7o 

(57) 

where  B,  is  the  ICL  bunching  parameter  and  6X  is  the  ICL  mismatch  factor.  The 
gain  formula  is  the  essencially  the  same  for  a  linear  polarized  ICL;  it  has  a  different 
coefficient  because  only  one  component  of  the  electrons’  transverse  motion  and  the 
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wave  field  (Eq.2*5)  are  applied.  B,  depends  on  the  beam  drift  velocity  while  Bc  is 
independent  of  it.  In  CARMs,  when  electrons  interact  with  the  wave,  some  electrons 
gain  energy  and  some  lose  energy  such  that  the  mismatch  factor  varies  in  time,  but  the 
bunching  parameter  always  remains  constant.  Thus,  the  electron  will  remain  at  the  same 
region  during  the  entire  interaction.  In  that  case,  change  of  the  mismatch  sign  may  cause 
the  saturation  of  wave  amplification.  However,  in  ICLs,  the  curve  dividing  regional) 
and  region(III)  depends  on  the  electron  axial  velocity.  During  ICL  interaction,  both  the 
mismatch  and  bunching  parameter  are  changing.  It  becomes  possible  for  electrons  from 
region(III)  to  enter  into  region(II)  or  vice  versa  after  they  have  lost  energy. 

From  Eq.  (41)  with  the  upper  sign  we  know  that  if  the  electron’s  vz  is  decreasing 
while  its  energy  is  decreasing,  the  resonance  interaction  can  be  prolonged  to  achieve 
higher  efficiency.  On  the  other  hand,  if  vx  is  increasing  while  its  energy  is  decreasing, 
there  will  be  no  autoresonance;  this  is  the  case  in  region(IV)  for  CARMs  and  ICLs. 
From  Eq.  (13),  we  note  that 

d  d 

=  ^Vg  ~  V°^dt7  (58) 

where  vg  =  c2/vph  is  the  group  velocity  of  the  wave.  The  first  term  on  the  right  hand 
side  of  Eq.  (58)  is  from  the  wave  magnetic  field  (the  v  x  B  force),  which  converts 
the  electron  axial  momentum  to  azimuthal  motion.  The  second  term  comes  from  the 
wave  electric  field  (the  mass  variation)  ,  which  increases  the  electron  axial  velocity  as 
it  dissipates  the  electron  azimuthal  energy.  While  the  electron  is  losing  its  energy,  and 
if  vg  <  voz(vg  >  voJ,  its  drifting  velocity  will  be  speeded  up  (slowed  down);  in  other 
words,  the  electron  always  increases  it’s  velocity  (in  absolute  terms)  relative  to  the 
wave.  Thus,  if  the  wave  group  velocity  (e.g.,  a  fast  wave)  is  slower  than  the  electron 
drifting  velocity,  that  is, 

c2  <  i'oIvph  ,  (59) 

then  the  change  in  axial  velocity  vt  (axial  bunching)  will  not  be  able  to  compensate  the 
change  in  energy  (azimuthal  bunching)  to  maintain  the  resonance.  The  competition  of 
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—* 

the  v  x  B  effect  and  the  mass  variation  effect  on  the  axial  velocity  also  exists  in  the 
bunching  parameter  and  the  rate  of  energy  exchange. 

In  a  circularly  polarized  FEL  configuration  (Fig.3c.),  the  helical  undulator  field  is 
represented  by 

Bw  =  By,(,x  sin  kwz  -  y  cos  kwz)  (60) 

where  Bw  is  the  amplitude  of  the  external  magnetic  wiggler  and  kw  is  the  wavenumber 
of  the  wiggler  field.  The  radial  driving  force  comes  from  the  vo*  x  Bw  term.  In  this 
case  Fo  =  'yomevo±kwvot  is  the  amplitude  of  the  driving  force,  u0x  =  (eBw)/(^omtkwc) 
is  the  electron  azimuthal  velocity,  q  =  0,  wt  -  kwc,  fo  =  kwzo  and  k  »  kw  for  energetic 
electron  beams  (v0t  — »  c).  The  resonance  condition  can  be  expressed  by 


u>  —  kvx  —  kwc  ~  0 


(61) 


which  means  that  the  resonant  wave  frequency  is  approximately  2^ltkwc.  The  rate  of 
electron  energy  change  and  the  gain  are 


e2E] t&w  j 

iw-'  B’6‘ 


„  .  IcosSft  -  2  +  6/tsin6y< 

Ga,n  ~  - S’ - 

Bf  =  0-(4-~ ) 
vi 


vB 


-p 

{/  =  w  -  kv Oz  —  ky,C 


(62) 

(63) 

(64) 

(65) 


where  Bf  is  the  FEL  bunching  parameter  and  6/  is  the  FEL  mismatch  factor.  The  gain 
result  can  be  recovered  from  Ref.  [18],  which  is  based  on  solving  a  pendulum  equation, 
by  letting  Bf  ~  For  an  FEL  (q  =  0),  the  gyrofrequency  of  the  electron  is 

independent  of  the  beam  energy  so  that  the  point  P  (i.e.,  uo,  — ►  0)  in  Fig.2.  is  at 
infinity.  Region  (III)  vanishes  and  there  is  no  azimuthal  bunching. 

Due  to  the  pondermotive  force  (Eq.  12),  the  interaction  of  the  wave  and  the  electron 
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in  the  axial  direction  is  changed  to 


i>o2.  d 


(66) 


where  vp  =  u/(k+kw).  We  can  rewrite  c(—  ~xu)  as  (ve+vep-v0z)  where  vgp  =  <?kwju}. 

Vp  C 

The  pondermotive  force  of  the  magnetic  wiggler  enhances  the  EM  wave  t*>  x  B\  force. 
For  practical  puiposes,  the  mismatch  is  made  very  small  so  that  the  electron  drifting 
motion  is  synchronous  with  the  pondermotive  potential  ( vp  ~  t’o2).  In  this  situation, 
{f-  —  *»)  ~  ^-(1  —  oc  and  is  always  positive.  Therefore,  the  change  of  the 
electron  axial  velocity  vz  has  the  same  sign  as  its  energy  change.  The  conventional  vxB 
force  is  always  larger  than  the  mass  variation  effect  in  the  axial  bunching  process. 

Table  II.  summarizes  the  comparison.  All  the  driving  forces  are  transverse.  Their 
resonance  conditions  are  similar  but  with  different  q  indices  on  the  gyrofrequencies. 
The  P  points  which  separate  region  (II)  and  region  (IV)  are  l,v2  and  oo  for  CARMs, 
ICLs  and  FELs,  respectively.  Their  rates  of  electron  energy  change  are  the  same  except 
for  different  bunching  parameter  and  mismatch  factor.  FELs  have  the  most  favorable 
scaling  for  the  output  frequency  with  the  beam  energy  (oc  72),  but  the  electron  undulation 
frequency  is  limited  by  the  available  wiggler  period.  Thus  it  requires  high  7  to  achieve 
high  frequency  lasing  at  the  expense  of  reducing  the  gain.  CARMs  have  the  best 
controls  on  beam  guiding  and  wave  mode  selection,  but  the  radiation  output  frequency 
is  only  linearly  proportional  to  7  and  the  electron  gyration  frequency  is  limited  by  the 
available  magnet  technology.  ICLs  can  have  the  highest  electron  oscillation  frequency 
(e.g.,  u>ei  ~  10* GHz  for  ~  10,6cm“3)  and  the  output  frequency  scales  more  favorably 
than  for  CARMs,  i.e.,  ~  7^.  In  addition  to  these  advantages,  ICLs  do  not  require  any 
applied  external  field. 


21 


IV.  Summary  and  Discussion 


We  have  developed  a  unified  small-signal  theory  to  describe  the  amplification  mech¬ 
anism  in  relativistic  electronics.  The  mechanisms  persist  to  nonrelativistic  electron  en¬ 
ergies.  This  study  unifies  the  amplification  mechanisms  of  different  wave  emission 
devices  by  assuming  that  the  electron  oscillation  frequencies  depend  on  their  beam  en¬ 
ergies  as  7-<J  where  q  =  1,|  and  0,  respectively  for  CARMs,  ICLs  and  FELs.  Their 
differences  basically  come  from  different  q  values. 

The  azimuthal  bunching  which  depends  on  the  q  value,  electron  drift  velocity  and 
wave  phase  velocity  is  induced  by  the  wave  electric  field  ( E i)  and  is  a  relativistic  mass 
effect.  On  the  other  hand  the  axial  bunching  caused  by  the  v  x  B  force  which  is  a 
nonrelativistic  effect  and  the  relativistic  mass  variation  is  independent  of  the  q  value. 
The  bunching  can  take  place  in  velocity  space  related  to  wave.  It  is  not  necessary  for 
electrons  to  be  bunched  in  real  axially  space  to  have  wave  amplification.  The  competing 
effects  which  cause  bunching  exist  in  both  the  single  electron  resonance  condition  and 
the  collective  rate  of  electron  energy  change.  The  azimuthal  bunching  is  dominant 
over  the  axial  bunching  if  q(l  -  >  (fr-^)»  and  vice  versa.  If  voz  >  vp  or 

c2  <  vozvp  ,  in  stead  of  offsetting  the  two  bunchings  reinforce  each  other  which  results 
in  no  autoresonance.  In  order  to  have  wave  amplification  the  sign  of  the  mismatch 
is  required  to  be  the  same  as  the  sign  of  the  bunching  parameter.  For  CARMs  and 
ICLs  the  relativistic  mass  variation  effect  on  the  electron  axial  velocity  may  overwhelm 
the  v  x  B  force  so  that  the  change  in  axial  velocity  is  opposite  to  the  electron  energy 
change  and  the  axial  bunching  term  becomes  negative.  For  FELs  the  bunching  arising 
from  the  relativistic  mass  variation  is  smaller  than  that  of  the  pondennotive  force  from 
the  wiggler  and  wave  magnetic  fields  so  that  the  FEL  bunching  parameter  is  always 
negative.  All  these  mechanism  can  upshift  the  undulation  frequency  by  a  factor  of  2-y^. 
The  ICL  electron  gyrofrequency  w,  =  (27rn,e2)/(me7l)  can  be  made  to  be  the  highest 
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among  the  three  with  the  present  day  technology.  FELs  and  CARMs  already  have 
successful  experimental  results.  ICLs  are  expected  to  become  a  realistic  mechanism 
to  generate  coherent  tunable  radiations  from  the  microwave  regime  to  the  soft  X-ray 
regime. 

The  unified  small-signal  theory  not  only  links  the  old  physics  issues  in  CARMs  and 
FELs,  but  also  reveals  the  new  physics  issues  in  ICLs.  Any  new  mechanisms  ought  to 
be  covered  by  the  unified  theory. 
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Table.I.  For  wave  emission  the  relation  between  bunching  parameter  and 
mismatch  factor. 


bunching  parameter  B 

required  mismatch  6 

dominant  bunching 

positive 

positive 

azimuthal  f 

0 

no  wave  emission 

cancel 

negative 

negative 

axial 

t  For  FEL  there  is  only  an  axial  bunching  and  no  azimuthal  bunching. 
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Figure  captions 


Fig.l.  A  basic  electron  relativistic  radiation  model. 

Fig.2.  Relations  between  two  bunching  mechanisms  in  four  vqz/c  vs.  vp/c  regions. 

Fig.3.  The  electron  motion  and  external  fields  of  CARM,  ICL,  and  FEL.  Their  external  forces 
play  a  similar  role  although  the  origin  and  direction  of  fields  are  quite  different. 
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